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Abstract 

We consider a random, uniformly elliptic coefficient field a(x) on 

the d-dimensional cubic lattice Z,'^. We are interested in the spatial 

Q> . decay of the quenched elliptic Green's function G{a; x, y). Next to sta- 

P^ I tionarity, we assume that the spatial correlation of the coefficient field 

r^ • decays sufficiently rapidly to the effect that a Logarithmic Sobolev 

Inequality holds for the ensemble (•). We prove that all stochastic 
moments of the first and second mixed derivatives of the Green's func- 
tion, that is, {\VxG{x,y)\^) and {\VxVyG{x,y)\P), have the same de- 
cay rates in |x — y| 3> 1 as for the constant coefficient Green's function, 
^ I respectively. This result relies on and substantially extends the one 

OO ■ by Delmotte and Deuschel [8], which optimally controls second mo- 

Zl^ ■ ments for the first derivatives and first moments of the second mixed 

T^ , derivatives of G, that is, {\VxG{x,y)\'^) and (|Va;VyG(x,y)|). As an 

t;-*^ I application, we derive optimal estimates on the random part of the 

^D ' homogenization error. 

cn 

The outline of this work is as follows: After introducing the discrete setting 
in Section 1, we present the statistical assumptions and the main result on 

K> ' the annealed moments of the Green's function in Section 2. We present opti- 

H ■ mal estimates on the random part of the homogenization error in Section 3. 

■ ■ ■ Section 4 contains an annealed Holder-estimate in the spirit of De Giorgi, 

which we will obtain as a a consequence of our main result. In Section 5 we 
explain our main assumption, a weakened Logarithmic Sobolev Inequality, 
which in particular holds for all independent, identically distributed coef- 
ficient fields. Section 6 contains the main ingredients of the proof of the 
annealed Green's function estimates — in particular we recall the result by 
Delmotte and Deuschel [8]. All proofs are postponed until Section 7. 
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1 Discrete uniformly elliptic equations 

In this paper we consider linear second order difference equations with uni- 
formly elliptic, bounded random coefficients of the form 

V*{aVu){x) = fix) forallxeZ'^. 

If there is no danger of confusion, we also write V*aVti for V*{aVu). In this 
equation we define the spatial derivatives for scalar fields ^ : Z'^ — )■ M, vector 
fields ,^ = (,^1 , . . . , ,^rf) : Z*^ — )■ M.'^, and alH = 1, . . . , d by the expressions 

V.C(x) := C(a; + ei) - C(x), V*C(a;) := ({x - e,) - C(x), 

d 
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Here ei, . . . , e^ is the canonical basis of M"'. The expressions VC and — V*^ 
are the discrete gradient and divergence for scalar fields and vector fields, 
respectively, on Z'^. Note that V* is the £^(Z'^)-adjoint of V. 

Here comes the main assumption on the considered coefficient fields {a{x)}^^id: 
For every x G Z"^, a{x) G R"^^"^ is symmetric and uniformly elliptic in the 
sense of A|^p < ^ ■ a{x)^ < |^p for all ^ G M'^. Here and below A G (0, 1) 
denotes the ellipticity ratio which is fixed throughout the paper. In addition 
we assume that the coefficients are diagonal. We note that the condition 
a = diag(aii, . . . , add) is diagonal yields that the coefficients can be associ- 
ated with the bonds of the lattice (rather than the vertices): In fact, this can 
be seen on the level of the Dirichlet integral, which can be rewritten as 

d 

Y, C(a:)V*(aVC)(x) = E E ^nix)iax + e,) - C(x))^. 

This allows, for instance, to interpret V*a\/u as the generator of a random 
walk on Z'^ with jump rates across bonds described by a. Diagonality is 
known to be a sufficient (but not necessary) condition for the maximum 
principle to hold for V*aVu. The maximum principle is a crucial ingredient 
for the estimate (19) on the quenched Green's function. Also, in the proof of 
Lemma 5, we shall adapt the bond-based point of view. 

Our main object is the non-constant coefficient, elliptic, discrete Green's 
function G{a;x,x') defined through V*a(x)VxG(a; a;, x') = 6{x — x'), where 
6 stands for the discrete version of the Dirac distribution, i.e. 

^. , I 1 for X = 
d[x) = < 

otherwise 



We usually drop the argument a and just write G{x,y). Often, it is more 
convenient to appeal to the distributional characterization: 

V C(x) : Yl ^^(^) ■ «(a^)V..G(x, x') = C(x'). 

Throughout the paper, we will work in dimension d > 2. Dimension d = 2 
needs a bit more care in terms of the definition of the Green's function. 
Since we are only interested in gradient estimates, this is merely technical. 
Sometimes, it is more convenient to think of the approximation via a massive 
term in the sense of 

T-^Gt{x, x') + VXa:)V^G'r(x, x') = 6{x - x'); 

this is the case in the proof of Proposition 1. At other times, it is more 
convenient to think in terms of an approximation via periodization in the 
sense of 

VXx) V^Gl(x, x') = Ys{x-x' - Lz) - L-'^- (1) 

this is the case in the proof of Lemma 6. 



2 Assumptions on the ensemble and main re- 
sult 

We are given a probability measure on the space of uniformly elliptic, di- 
agonal coefficient fields (endowed with the product topology), cf. the previ- 
ous section. Following the convention in statistical mechanics, we call this 
probability measure an ensemble and denote the associated ensemble av- 
erage (i.e. the expected value) by (■). We assume that (■) is stationary 
in the sense that for any shift vector z G Z'^, the shifted coefficient field 
a(- -|- z) := {U^ 3 X ^^ a{x + z)} has the same distribution as a. We also 
note that the Green's function is shift-invariant or stationary in the sense 
that G{a{- + z);x,y) = G{a] x + z,y + z). 

Besides stationarity, the main assumption on the ensemble of coefficients 
and only probabilistic tool will be a weak variant of the Logarithmic Sobolev 
Inequality (LSI). In Section 5, we will comment on the LSI and the related 
Spectral Gap Inequality — there we will also describe the relation between 
this weak LSI and the usual LSI. 

Definition 1. [Weak logarithmic Sobolev inequality]. Let (■) be a stationary 
ensemble of coefficients a. 



i = 1,- ■ ■ ,d, we define the continuum vertical derivative g^j . as the partial 
derivative with respect to the variable aii{y). Furthermore, we denote by q^S^ 
the d-dimensional vector with coefficients „ % . with i = 1, . . . ,d. Thus in 



part^cular,i^y = EU^S^r■ 

a) We say (■) satisfies a weak Logarithmic Sobolev Inequality (LSI) with 

constant p > if for all positive random variables f{a) > we have 

Note that the difference between the wealcened LSI (2) and the usual LSI, see 
(15), lies solely in the supremum over a{y). The usual LSI therefore implies 
(2). The merit of this weakening is that it is satisfied by any ensemble 
of independent, identically distributed coefficients {a{y)}y^^d, cf. Lemma 1 
below. Our main result is: 

Theorem 1. Let (■) be stationary and satisfy the weak version of LSI with 
constant p > 0, see Definition 1. Then for all 1 < p < oo and x G Z'^, it 
holds 

(|VVG(x, 0)1^0^ < CKA,p,p)(|a;| + l)-^ (3) 

(|V.G(x,0)p^)^ < C{d,X,p,p){\x\ + lY-'. (4) 

Here VVG(x, y) = Vx^yG{x, y) denote the mixed second derivatives of the 
Green's function. Also, here and in the sequel, C{d, X, p,p) stands for a 
generic constant that only depends on dimension d > 2, on the ellipticity 
ratio \ > 0, on the LSI constant p > 0, and on the exponent of integrability 
p < oo. 

Clearly, the spatial decay rates in Theorem 1 are optimal, since those are the 
decay rates of the constant coefficient Green's function. Note that station- 
arity of (■) and G implies (|VVG(a; x,|/)|2p) = (|VVG(a(- - y); x,i/)pP) = 
(|VVG(a; x — y, 0)^^), so Theorem 1 can be rephrased as 

(|V,V,G(x,2/)p^)i < C(rf,A,p,p)(|x-y| + l)-^ 
(|V.G(x,2/)p0i < Cid,X,p,p){\x-y\ + iy~'' 

for all x,y & JJ^. An interesting aspect of Theorem 1 is the following: The 
quenched versions of (3) and (4) are false, i.e. the uniform in a and point- 
wise in X estimates |VVG(x,0)| < G{d,X){\x\ + 1)"'^ and |V^G(x,0)| < 




C((i, A)(|x| + I)'' ^ do not hold (while suitably spatially averaged versions of 
both estimates do hold uniformly in a); see our discussion in Section 4 below. 

An easy consequence is the following generalized variance estimate on G 
itself: 

Corollary 1. Let (■) be as in Theorem 1. Then it holds that 

I 9r,\ - fflrl + 112(1-'^) 

\ / l(|x| + l) ^log(|a;|+2) d = 2 

jor all X ^ Tj'^ and 1 < p < oo. 

Remark 1. We note that the estimate in Corollary 1 is optimal in the scaling 
of the spatial decay. This can he seen by developing to leading order in a small 
ellipticity ratio 1 — A ^ 1. We expand upon this argument (for the special 
case of p = 1) after the proof of Corollary 1. 



3 Homogenization error 

In the same vein as Corollary 1, Theorem 1 allows to give optimal estimates 
on the random part of the homogenization error. These extend the results 
by Conlon & Naddaf [5, Theorem 1.2, Theorem 1.3] from small ellipticity 
ratio (i.e. 1 — A -C 1) to arbitrary ellipticity ratio. For the "strong error" (see 
below for an explanation of this wording) [5, Theorem 1.2] in (i > 3, this was 
already achieved by Gloria [11, Theorem 2]. For all other cases, our result 
appears to be new. Let us be more precise: For a coefficient field a{x) and a 
right hand side f{x) we consider the solution u{x) of 

V*a{x)Vu{x) = fix) on Z'^. (6) 

In order for (6) to have a unique solution that decays (i. e. lim|a.|i-ooM(x) = 0), 
we assume for simplicity that / is compactly supported (and is of zero spatial 
average in the case of rf = 2 to ensure decay of u). By the random part of 
the homogenization error, we understand the "fluctuations" u{x) — {u{x)). 
These are expected to be small (w. r. t. the size of u{x) itself) if f{x) varies 
only slowly w. r. t. to the lattice spacing. In our notation, the lattice spacing 
is unity, so that a natural model for a right hand side that has a large 
characteristic scale L S> 1 is given by /(x) = L~^/(^) for some bounded and 
compactly supported "mask" f{x), x G M.'^. The scaling L~^ of the amplitude 



of / is motivated as follows: In the rescaled variables x, (6) now assumes the 
suggestive form of 

VXf )V,m(x) = /(x) oneZ^ (7) 

where e := L~^ is the ratio of the lattice spacing to the characteristic scale 
of the r. h. s. and where V^ denote finite differences for the rescaled lattice 

eZ'^ (i. e. Ve,jM(x) = e^^{u{x + ecj) — u{x))). 

The size of the fluctuations will be measured in two different ways. 

• Corollary 2: Here, the fluctuations will be controlled in a strong way in 

the sense that we estimate the (discrete, spatial) £P(Z"')-norm {J2x I"" ~ ('")D 
of the fluctuations. This will be done for arbitrary stochastic moments 
(the role played by rp). Corollary 2 is the generalization of [5, The- 
orem 1.2] as well as [11, Theorem 2]. For our model right hand side, 
f{x) = e^/(ex) with bounded and compactly supported /, the fluctua- 
tions are (up to a logarithmic correction for d = 2) oi the order of e in 
this measure, see (11). 

• Corollary 3: Here, the fluctuations will be controlled in a weak way 
in the sense that we only estimate spatial averages ^^.(m — {u))g of 
the fluctuations, with deterministic averaging function g{x). Again, 
this will be done for arbitrary stochastic moments (the role played 
by r). Corollary 3 is the generalization of [5, Theorem 1.3]. For our 
model right hand side f{x) = e^f{ex) with bounded and compactly 
supported /, and an averaging function of the form g{x) = g{ex) with 
bounded and compactly supported g, the fluctuations are 0(e'^/^) in 
this measure, see (12). 

Corollary 2. Let (■) be as in Theorem 1; for compactly supported right 
hand side f{x) consider the decaying solution u{x) to (6). Let the spatial 
integrability exponents 2 < p < oo and 1 < q < oo be related through - = 

In case of d > 2, we have for all r < oo.' 

(Eh-(")r)y'^C'(d,A,p,p,r)(5^|/|' 

^ X ^ ' X 

In case of d = 2, we additionally require p > 2 (so that q > 1) and that f is 
supported in {\x\ < R} for some R> 1. Then we have for all r < oo.' 



\x\<R 




El^-(^)r ) <C(A,P,P,r)(log^i?)(5^|/r)l (9) 



Corollary 3. Let (■), f{x), u{x) be as in Corollary 2. Let the averaging 
function g{x) be compactly supported. Let the two integrability exponents 
1 < q,q < oo be related ^?/- + i = | + |- Then it holds for all r < oo: 



Y.(^-{^))9 y <C{d,X,p,r)lj2\fn(J2\9 



(10) 



For the convenience of the reader we express the results of both corollaries 
in terms of the rescaled variable x = ex, the model right hand side f{x) = 
e^/(ex) and the model averaging function g{x) = g{ex); we also rewrite the 
solution itself in terms of u{x) = u^iex). In this notation, (8) (multiplied by 
^d/p-j \^■^J^Yla.s into 



< C(rf,A,p,p,r)e(e^5^ |/>]' < C(rf, A,p,r, />. (11) 

Note that this can be interpreted as the discrete version of 

lu^-iuM^dx) ) <C{d,X,p,p,r)e{ I Ifl^dx] 



which highlights the 0(e)-nature of the "spatially strong" error. 
Likewise, (10) (multiplied by e"^) turns into 



e'Y.{u,-{u,))g 



xGel 



< C{d,\,p,r)eUe''Yl 1^1 



xGel 



in le'^J2 \9\' 



xGel 



< C{d,X,p,rJ,g)e2. 
As above, this can be seen as the discrete version of 



{u, - {u,))g 



dx) < C{d, A,p, r)e^ 



\f\'dx 



(12) 



\g\^dx 



uncovering the 0(e'^/^)-nature of the "spatially weak" error. 



Let us make a couple of further more detailed remarks related to Corollaries 
2 and 3. In case of Corollary 2 and d = 2, we can use Holder's inequality to 
establish an estimate also for p = 2. However, in that case we pay the price 
of an arbitrarily small power of R on the right hand side of (9). We also 
note that the requirement that / has compact support and that u decays 
can be weakened: All we need is the Green's function representation u{x) = 
J2y^{^^y)f{y)- We conclude by pointing out that our argument does not 
require any smoothness assumptions on f[x) and g{x) beyond (uniform) 
boundedness to obtain (11) and (12). 

In its e-version (12), the estimate of the weak error (10) hints towards the 
scaling of the Central Limit Theorem (CLT). CLT-scaling has been estab- 
lished for related quantities, as we shall detail now: For g = f, the weak 
measure of fluctuations turns into a measure of fluctuations of the energy: 

y ^{u — (u))g = 2^ Vn ■ aVu — (y^ Vu ■ aVu). 

XXX 

If u is the so-called corrector (which is an a-harmonic function of affine behav- 
ior on large scales) the (stationary) energy density deflnes the homogenized 
coefficient. In [12, Theorem 2.1], it is shown that in the case of indepen- 
dent, identically distributed (i. i. d.) coefficients, the energy density of the 
corrector has CLT scaling in the sense that spatial averages behave as if the 
energy density was independent from site to site; in [13, Proposition 7], that 
has result has been generalized to ensembles that only satisfy a Spectral Gap 
condition. The scaling result has been substantially sharpened for i. i. d. 
ensembles: In this situation, the fluctuations of the energy of the corrector 
become more and more Gaussian as the box over which the spatial average 
is taken increases. The latter result has been obtained by three different 
techniques: Nolen [23] gives a quantitative estimate based on a differential 
characterization of Gaussian distributions (second order Poincare inequality) 
and relies on the corrector estimates from [12, Theorem 2.1]. Biskup, Salvi & 
Wolff [3] obtain a more qualitative result using a Martingale decomposition 
of the spatially averaged energy density (their result assumes small ellipticity 
contrast 1 — A ^ 1, but presumably could be extended using the results of 
[13]). Rossignol [25] in turn uses an orthogonal decomposition of the space 
of coefficients (Walsh decomposition). 



4 Relation to De Giorgi's approach to elliptic 
regularity 

While our result heavily relies on the celebrated regularity theory for scalar 
elliptic operators, connected with the names of De Giorgi, Nash, and Moser, 
it also gives a new perspective on these results. We will specify the input from 
regularity theory, namely Nash's (upper) bounds on the parabolic Green's 
function, in the next section. We now address what we see as a new perspec- 
tive on these results, namely on De Giorgi's result on Holder continuity of 
a-harmonic functions. 

An elementary consequence of the mean value property is the following Liou- 
ville principle: Harmonic functions that grow sub-linearly must be constant. 
This holds for the constant-coefficient Laplacian both on Mf^ and on Z*^, but 
is no longer true for variable coefficients, even if they are uniformly ellip- 
tic. Indeed, a well-known example [1, Corollary 16.1.5] shows that for any 
a > 0, there exists an explicit coefficient field a^ < a{z) < 1 such that 
u{z) = TZe{\z\°'~^ z) is a-harmonic in z 3 C = M^. We believe that this exam- 
ple can be adapted to the lattice Z^ (provided the condition of diagonality 
is relaxed to the condition that the discrete maximum principle is valid, a 
setting to which our results presumably can be extended). However, a cele- 
brated result of De Giorgi [6, Theorem 2] states that for any dimension d and 
any ellipticity ratio A, there exists an exponent ao{d, A) > with the follow- 
ing property: For any field of coefficients X < a{x) < 1 and any a-harmonic 
function u{x), \u{x)\ < jxl"" for |x| ^ 1 implies that u is constant. This 
result holds both in R'' and in Z'^ [7, Proposition 6.2]. In this sense, while 
it is no longer true that "sub-linear implies constant", it remains true that 
"very sub-linear implies constant" . 

De Giorgi's result can be rephrased as an inner regularity result in terms of 
Holder continuity with Holder exponent Oq: For any harmonic function u{x) 
on {x : |x| < R}, the Holder-ao modulus of continuity at zero is estimated 
by the supremum: 



sup ^^^^^^-^^<C(rf,A)i?-"° sup \uix)\. 

x:\x\<R \Xr° x:\x\<R 

To contrast De Giorgi's result with our result below, let us rephrase it as 
follows: 

\u{x)-uiO)\ 

VA<a(x) < 1, Vi?<oo: sup , "'"'-"^ 7^^<C(rf,A), (13) 



where the outer supremum is taken over all u{x) that satisfy V*aVu = in 
{x : |x| < R}. 

In this context, we will show that Theorem 1 has the following Corollary. 
Corollary 4. For all < a < 1, p < oo, and R < oo, we have 

^^P- Tf^ ] ) < C'(^'^'"'P)' (14) 

u -^sup3,.|3,|<^ |-u(^a;j|y / 

where the outer supremum is taken over all u{x) that satisfy V*aVM = in 
{x : |x| < R}. 

As a consequence of this annealed Holder regularity, we obtain the following 
Liouville principle. 

Corollary 5. For {■) -almost every a, if u is a solution to V*aVu = which 
is strictly sub-linear in the sense that there exists < a < 1 such that 

lim -B^ sup \u{x)\ =0, 
then u is constant. 

Loosely speaking. Corollary 4 implies that for "most" coefficient fields, an 
a-harmonic function u{x) is Holder continuous with an exponent arbitrarily 
close to one. More precisely, the modulus of near-Lipschitz continuity of 
u{x) in some large ball {x : \x\ < R} only depends on its supremum in the 
ball of double radius. This quantitative result has the Liouville principle as 
an easy corollary: For almost every a, any sub-linear a-harmonic function 
must be constant, see Corollary 5 for the version of this result we can derive 
from Corollary 4. For the convenience of the reader, we display the short 
proof of how Corollary 4 yields Corollary 5. However, surprisingly for us, 
the qualitative Corollary 5 holds without any assumption on the ensemble (■) 
besides stationarity! This is established in a very inspiring paper [2, Theorem 
3]. The main ingredients for the short and elegant argument are 

• The "annealed" estimate (^^ |xpG(t, x, 0)) ^ t on the second moments 

of the parabohc Green's function G(a; t, x, y) = G(t, x, y) (cf. [2, 
(SBD)], see Subsection 6 below for the definition of G), which in our 
uniformly elliptic context even holds in its stronger "quenched" version, 

that is, E^|a;pG(t,a;,0) <t. 

10 



The annealed estimate —{J2x ^i^^ ^' 0) liiG'(t, x, 0)) < logt on the spa- 
tial entropy of the parabolic Green's function G (cf. [2, p. 12]), which 
in our context is an immediate consequence of the second moments 
estimate. This ingredient is shown to imply the following annealed 
continuity property of G: 



^^ ' '^^^ Git,0,x) + Git-l,y,x 



x) I ^ t 



5 Logarithmic Sobolev inequality 

In the following we give a more detailed description of our use of the log- 
arithmic Sobolev inequality and prove that any i. i. d. ensemble satisfies 
Definition 1. LSI substitutes the Spectral Gap Inequality (SG) in prior work 
on quantitative stochastic homogenization. SG has been introduced into the 
field by Naddaf & Spencer [18, Theorem 1] (in form of the Brascamp-Lieb 
inequality) and used most recently in [12, Lemma 2.3] in an indirect way and 
in [13] explicitly. Like SG, LSI follows from the property that there is an 
integrable fall-off of correlations in the sense of a uniform mixing condition 
a la Dobrushin-Shlosman, see for instance [26, Theorem 1.8 c)] for a discrete 
setting. Both SG and LSI quantify ergodicity of the ensemble, see for in- 
stance the discussion in [13, Chapter 4]. Recall that the usual LSI in this 
setting (with continuum derivative) would read 

The difference to Definition 1 lies solely in the presence of a supremum over 
a(y) inside the expectation, see (2). 

Both SG and LSI are based on the notion of a vertical derivative that defines 
a Dirichlet form and thus a reversible dynamics, namely Glauber dynamics, 
on the space of coefficient fields (the word "vertical" is used to distinguish 
this derivative from the "horizontal" derivative naturally arising in stochastic 
homogenization, but not used in this paper). In the earlier work on stochastic 
homogenization and motivated by field theories, see [19], the version of SG 
that is based on the continuum vertical derivative (as on the r.h.s. of (15)) has 
been used [18]. However, this assumption rules out the natural example of 
coefficients with a single-site distribution that only assumes a finite number 
of values (Bernoulli). In [13], it has been shown that SG based on a discrete 
version of vertical derivative can be used instead — without much more work 
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and allowing for any single-site distribution. Indeed, it is known that SG for 
the discrete vertical derivative holds for any independently and identically 
distributed field of coefficients: This follows from the tensorization principle 
and the fact that any single-site distribution satisfies SG for the discrete 
vertical derivative with constant one. 

Here lies a difference between SG and LSI: While the tensorization principle 
still holds for LSI, see for instance [15, Theorem 4.4], it is not true that the 
appropriate discrete version of the usual LSI holds for any single-site distri- 
bution — in fact, the discrete LSI does not hold for a continuum distribution. 
Hence in order to treat arbitrary single-site distributions, we are forced to 
consider the weakened version of LSI of Definition 1. A similarly weakened 
version of SG was already considered in [12, Lemma 2.3]. 

The LSI has been of great use in the setting of stochastic processes and 
diffusion semi-groups, for the first time introduced in generality by Gross [14]. 
It implies SG and is equivalent to the notion of hyper-contractivity, see [14, 
Theorem 1] as well as [15, Theorem 4.1] for a recent exposition. Incidentally, 
hyper-contractivity was first observed in the Gaussian context by Nelson 
[21], see [22] for an improved result. It is thus the older notion and in fact 
motivated the (somewhat implicit) introduction of LSI by Federbush [10]. 
We refer to [15] for a recent exposition on LSI. 

The result of this section is that any independent, identically distributed 
coefficient- field satisfies the LSI (2) in Definition 1. 

Lemma 1. Consider an ensemble (■) of i. i. d. coefficients with arbitrary 
single-site distribution on [A, 1]*^. Then (2) holds, i.e. 

(/)/ ^p\^ a{y)e[x^Y f da{y)' I 

for all (continuously differentiable) positive functions f of the coefficient field 
a. The constant p can be taken to be p = ■^. 

Note that even for discrete distributions, we have to take the supremum over 
the whole box [A, l]'^, even though the coefficients attain only a countable 
number of values. Lemma 1 is an immediate consequence of the following 
two lemmas. The first one shows that any single-site distribution on [A, 1]*^ 
satisfies the LSI in Definition 1. 

Lemma 2. Let (■) be any distribution on [A, 1]'^. Then it holds that 

X\<1 sup J-(^ 
{f)' ~ 2p ae[A,i]d /(a) da 



/log^)<i sup ^(^r (16) 
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for all functions f : [A, 1]"' — )■ (0, oo). In fact, the constant p = ^ will do. 

The next lemma shows that the LSI in Definition 1 satisfies the tensorization 
principle. 

Lemma 3. Let (■) be an ensemble consisting of independent single-site dis- 
tributions such that each single-site distribution satisfies the single-site LSI 
(16) with the same constant p. Then {■) itself satisfies the LSI (2) with 
constant p. 

6 Main ingredients of the proof 

Loosely speaking, our approach consists in upgrading the (optimal) annealed 
estimates of Delmotte & Deuschel [8, Theorem 1.1] in terms of the integra- 
bility p: 

Proposition 1. [Delmotte & Deuschel] . Let {■) be stationary. Then we have 
for all X ^U^ 

(|VVG(x,0)|) < C{d,X){\x\ + l)-\ (17) 

(|V,.G(x,0)|) < C{d,\)i\x\ + lY~''. (18) 

More precisely, we refer to the estimates (1.4) and (1.5a) in [8, Theorem 
1.1] on the discrete parabolic Green's function G(t,x,y) = G{a;t,x,y) (i.e. 
the solution of dtG(t,x,y) + Vla{x)'VxG(t,x,y) = with G(t = 0,x,y) = 
6{x — y)) that in our notation imply for any weight exponent a < oo: 

(|VVG(t,x,0)|) < C(rf,A,a)(t + l)"^i(^ + l)"^ (19) 

(|V.G(t,x,0)|) < C(rf,A,a)(t+l)-f-^(J^ + l)"l (20) 

(In fact, [8] establishes (19) & (20) with exponentially decaying weights in- 
stead of just algebraically decaying ones.) Since the elliptic Green's function 
can be inferred from the parabolic one via G{x, y) = J^ G{t, x, y)dt, these es- 
timates imply (17) & (18) (by fixing some a > d and performing the change 
of variables t = \x\~'^{t -\- 1)). Actually, [8] establishes (20) and thus (18) 
in the stronger form where the L^-norm (| ■ |) is replaced by the L^-norm 
(I ■ \y/^: (I V.G(x, 0)p) V2 < c{d, A, a){\x\ + l)i"'^(g; + l)"f . 
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Let us point out that the spatially point-wise annealed estimates (19) & (20) 
are consequences of the following spatially averaged quenched estimates 

J2{{^ + l)^G{t,x,0))' < C{d,\,a){t + 1)-I, (21) 

X 

E((^ + 1)^I^-^(*'^'0)I)' ^ Cid,\,a)it+irl-\ (22) 

X 

The first estimate (21) is the (upper, off-diagonal part of the) celebrated Nash 
estimate [20, Appendix]. The discrete case was treated in full generality in 
[4, Corollary 3.28]. The second estimate (22) is a consequence of the first 
one. For an elementary proof of both, we refer to [13, Lemmas 24 & 25], with 
the Nash inequality as only noteworthy ingredient. Let us point out how (22) 
is implies (19): Using the semi group property in form of WG{t,x,x') = 
J2y^xG{^,x,y)'Vx'G{l,y,x') we obtain by the triangle inequality for the 
weight, Cauchy Schwarz in ^ , and the symmetry of G{t, x, y) in x and y: 

^ E (^ + 1) V.G(i,.,.)i(i^ H- i)'iv.G(|.,.',i 



y 

xE((^Trf^ + i)*iv.-G(5.<j')i)7 



Note that the right-hand side of the last inequality does not allow for appli- 
cation of (22), since the sum is not in the variable the derivative is taken. 
However, taking the expectation, using Cauchy Schwarz in (■), and using sta- 
tionarity and symmetry in form of (I Va;G'(|,a;, y)p) = (|V2G(|,a; — y, 0)p) = 
(|V3G(|,0,x-?/)p) = (|VyG(|, -X, -y)p) (where V2 stands for the deriva- 
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tive w. r. t. to the second variable), yields 

l)^(|VVG(t,x,x')|) 



I /I2 

\x — x\ 



t+1 



y 



''l^ + l)"(|V..G(i 



xE(^^TTr- + ini^"''(^'"''!''l'> 



y 



^ y 



X 

y 



We now see that (22) implies (19). The estimate (20) is derived via the semi 
group property in form of Va;G(t, x,x') = ^ Va;G(|,x,|/)G(|,|/, x') from 
the combination of (21) and (22) by an analogous argument. 

Note that the estimates of Proposition 1 make no assumptions on the en- 
semble besides stationarity. In order to pass from Proposition 1 to Theorem 
1, we need the assumption on the ensemble from Definition 1. In fact, LSI 
enters only through the following lemma. 

Lemma 4. Let (■) he stationary and satisfy LSI with constant p > 0. Then 
for arbitrary 6 > and 1 < p < oo and for any ({a) it holds that 



a{x) 



In order to make use of Lemma 4, we need to estimate the vertical derivatives 
of VVG and VxG. The following lemma is at the core of our result. 

Lemma 5. There exists an integrability exponent po = po{d, X) < oo such 
that for all p > Po and any x G Z'^, it holds 

(|a:| + l)V(5^sup||-VVG(x,0)|T)'^ 

< Cid, A,p) sup jd^l + 1)^(|VVG(^, 0)pP)^} (24) 
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and 

<C(ci,A,p)(sup{(|z| + l)'^-i(|V.G(z, 0)1^0^} (25) 

+ sup|(|2| + l)^(|VVG(2,0)|'O^})- (26) 

The formulation of Lemma 5 shows that with our method, we first have to 
estimate the mixed second derivatives (|VVG(x, 0)^^) before we can tackle 
the first derivatives (|V3;G(a:, 0)^^). It also reveals that it is necessary to 
estimate high moments p > Pq in (■) in order to estimate moderately low 
moments like the fourth moment (|VxG(a;, 0)|^) that is needed in the proof 
of Corollary 1. 

The preceding lemma relies on the following suboptimal, but quenched esti- 
mates on the (elliptic) Green's function: 

Lemma 6. [Gloria & Otto] There exists «o = o^oid, A) > such that for all 
R>0 

|VVG(0,0)|2 + i?2"o J2 |VVG(x,0)|2 < C{d,X), (27) 

i?,<|x|<2_R 

|V,.G(0,0)|2+ Yl |V.G(x,0)p < Cid,\)- (28) 

R<\x\<2R 

The estimate (28) was established in the stronger (dimensionally optimal) 
form of '^ji<M<2R l^xG'(x, 0)p < _R^~'^ in [12, Lemma 2.9]; in its weaker form 
of (28), it is straight forward for d > 2. The proof of estimate (28) in [12] in 
case of d = 2 is subtle and relied on an adaptation of [9]. In this paper, we will 
give an elementary argument for the estimate (27), which we could not find 
in the literature. However, in the stronger form of J2r<\x\<2R |VVG'(x, 0)p < 
C{d, A)-R^~'^~^"°, this estimate can also be seen as a consequence of the fol- 
lowing classical ingredients: 

• the optimal decay of G{x, y) itself, that is just needed in a spatially av- 
eraged sense of R~^ ^yR<\x-v\<2R l^(^' v) ~ G\ < C{d, X)R^~'^ (thanks 
to subtracting the average G over the annulus {y : R < \x — y\ < 2R}, 
this estimate also holds in c? = 2), 



De Giorgi's Holder continuity estimate, that then yields for some Oq 



ao{d,\) > that sup,^^<u_,|<2/j |V,G(x,2/)| < C(d, A)/?^-^--^, 
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• Caccioppoli's estimate, that then yields J2x-R<\x-y\<2R \^y^xG{x, y) p < 
C(d,A)i?2-<i-2"o. 

Remark 2. We remark that with the same proof, one obtains a periodic 
version of Theorem 1 (with constants uniform in L) for the Green's function 
defined in (1). In that case, one just replaces the Euclidean distance \x\ 
on if- by its periodic version dist(x, LZ*^) on the torus M/LZ'^. The periodic 
version of Proposition 1 follows as above from the quenched spatially averaged 
estimates of [13, Theorem 3(b)]. 

7 Proofs 

Proof of Lemma 4. 

Step 1. Result for p = 1. We claim that for any 5 > and all C(o), it holds 

where p denote the constant in the LSI, see Definition 1. By homogeneity, 
we may assume (C^) = L For all real- valued C, it holds that 



C^< 



rexp(J,)|C| if|C|<expJ, 
\£fC'logC^ if|C|>expJ, 



Since x log x is bounded from below by -, it is elementary to verify that - |C| + 
C^logC^ > whenever \C^\ < 1. Since furthermore C^logC^ > whenever 
\(\ > 1, we find that 

Hence taking the expectation (■) yields 



\p6'^J 2eJ''^" 4 
Since (C^) = 1, Young's inequahty yields 



-''l^j^^j<i^|>'^+H^="'i^^i^i <^'>- 
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Combining the last two estimates, we deduce 



«^)^(-(^)^a<ici)^^^(c'- 



Hence LSI yields 

/ 2 \ p5 



(OS|exp^-j 



.g<,c,)v.XE-(^) 



and estimate (29) follows from taking the square root and applying the esti- 
mate A/a+~fe < y/a + y/h. 

Step 2. We finish the proof of (23), i.e. show that 

(C->*<C(p,M)(|CI).<((E-P(5|))T))* 

for general p > 1. To that end, we apply (29) to (^ replaced by |C|^ and take 
the 2p-th root: 



a[x) 



"da{x) 



where C(p,p,S) denotes a generic constant only depending on p, p, and 
5. Since p < 2p, an application of Holder's inequality in (■) and Young's 
inequality on the first r. h. s. term yields 

(|CrO^<C(p,p,5)(|C|) + 2(<5(^sup(^|Cr)'))'^. (30) 



X '^(^) 



Now we apply the chain rule q3^\C\^ = p\C\^ ^C^Ji^yC) since p > 1. Fur- 
thermore we note that for every coefficient field a E [-^j 1]^ , it holds that 
suPa(^) Id < ICI + Wsup„(^) \g^\ and hence 

!;2?(|^|*ia|))V'^(''-)(i^|'^'-^!r.'?(5|))'^(!r5(a^)T)' 

Hence 

(i:-p(5|;)Klf) 

.CK,((E,C,-=;P(^)V(E(-(^)T 

^-<^-'(<l^l^'>'-K(E:^p(a|y)T)^-((?-(^^^^' 
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Here in the last line, we have used Holder's inequality in (■) with exponents 
(2^,p) on the first term and the embedding £^(Z'^) C £^^(Z"') on the second 
term. Another application of Young's inequality thus allows us to pass from 
(30) to 

(|Cp^)^<C(rf,p,p,5)(|C|) + C(p)((5 + 5^)((5^sup(^)y))^. 

By redefining 6, this yields (23). 

Proof of Lemma 6. 

We just give the proof of (27); for (28), we refer to [12, Lemma 2.9]. 

Step 1. In this step, we derive the a priori estimate 

5^|VVG(x,0)p<C(rf,A). (31) 

X 

Indeed, note that the weak formulation of the defining equation for G reads 
V ax) : Yl ^^(^) ■ «(^) VxG(a:, x') = ({x'). 

X 

Taking the derivative w. r. t. the variable x' in some direction i = 1, ■ ■ ■ ,d 
yields 

V C(x) : Yl ^^(^) ■ «(a;) V. V.,,.G(x, x') = V.C(x'). (32) 

X 

The choice of ({x) = Vi^x'G{,x,x') (we address the question of admissibility 
of this test function below) yields 

^ Va;Vi,^.'G'(x,x') ■ a(a;) V^.Vi,^.'G'(a;, x') = Vi,a;Vi,^'G'(a;,a;')U=x'- 

X 

Since a{x) > A, this implies (31) in the explicit form of 

$^|V.V.,.,G'(x,x')P<A-l (33) 



We now turn to the question of admissibility of ({x) = Vi^x'G{x,x') as a 
test function for (32), i.e. the question of decay as |x| f oo of this function 
and its gradient. This issue can be circumvented as in Step 3 below through 
approximation by the periodic problem. More precisely, we consider the 
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periodic discrete elliptic Green's function Gl{x,x') = GL{a,x,x') of period 
L. Up to additive constants, it is characterized by the weak equation 

5^ VC(x) ■ a(x) V,.Gi(x, x') = C(x') - L"'^ Yl ^(^ 






xe[-§,^r -6[-*-#)'' 



for all periodic C(^)- With the same argument as above, we obtain 

5^ |V.V.,..Gl(x,x')P<A-2. (34) 

Since Gl(x,x') converges point-wise to G{x,x'), the latter implies (33) in 
the limit L 'j' oo by Fatou's lemma. Incidentally, lim^i-oo VxGl{x, x') may be 
taken as a definition of VxG{x, x') in the case of (i = 2, where G itself is not 
unambiguously defined. 

In the following steps, we use the fact that u{x) = Vi^x'G{x,x')\x'=o, i = 
1, ■ ■ • ,d, is a-harmonic away from x = to show that there exists a decay 
exponent ao{d, A) > such that for all R > C{d) we have 

J2 |VM|'<C(rf,A)i?-2°° J2 \Vu\^. (35) 

x:\x\>R a;:|a;|>l 

Together with (31), this implies (27). In Step 2, we'll formally treat the con- 
tinuum whole-space case. In Step 3, we will show how to make the continuum 
case rigorous by approximation through the continuum periodic case. More 
precisely, using (31), we will directly prove the estimate (27) in form of 

J2 |VM|'<C(rf,A)i2-2°o. (36) 

x:\x\>R 

In Step 4, we indicate the changes necessary to treat the discrete case. 
Step 2. Formal derivation of the continuum version of (35), that is 



f \Vu\'^dx<C{d,X)R-^''" f \Vu\^dx 

J{\x\>R} >''{|a.i>l} 



(37) 



for R> 1 and a function u{x) satisfying 

-V^ ■ a(x)V:,M = in{|x|>l}. (38) 

Indeed, let 77 (x) be a cut-off function for {x : |x| > 2R} in {x : |x| > R}. 
We test (38) with C, = rf{u — u), where u is the spatial average of u on the 
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annulus {x : R < \x\ < 2R}. It is a priori not clear that this is an admissible 
test function for (38); we shall address this in the next step. We appeal to 
the identity 

V(?7^(m — u)) ■ aWu = V{ri{u — u)) ■ aW{ri{u — u)) — (m — uYVrf ■ aVr], (39) 

which in view of A < a{x) < 1 turns into the inequality 

V(r/2(n - u)) ■ aVu > A|V(r/(n - u))\^ - {u - n)^|Vr/|2. (40) 

Hence from testing (38) we obtain 

A / |V(r/(M - u))\'^dx < {u- uf\Vri\'^dx, 

which by the choice of rj yields the Caccioppoli estimate 

[ \Vu\^dx<C{d,\)R~^ f {u-ufdx. (41) 

J{x:\x\>2R} J {x:R<\x\<2R} 

By Poincare's estimate on {x : R < \x\ < 2R} with mean value zero, this 
turns into 

/ \Vu\'^dx<C{d,X) I \Vu\'^dx, 

J{x:\x\>2R} J{x:R<\x\<2R} 

which can be reformulated as 

f \Vu\^dx<C{d,\) f \Vu\^dx, (42) 

J{x:\x\>R} J{x:R<\x\<2R} 

where C{d, A) as always denotes a generic constant only depending on d and 
on A whose value may change from line to line. 

A standard iteration argument now leads from (42) to (37): Introducing the 
notation 4 := /|^.|^|>2fc} iVupda:, (42) reads 

VA;g{0,1,---} 4<C(rf,A)(4-4+i), 

which with help oi 6 = 6{d, X) := 1 — ^, < 1 can be reformulated 

VA;g{0,1,---} 4+i<^4, 

or with help of Oq = ctoid, A) := ^^^ as 

VA;g{0,1,---} 4 < ^^/o = (2^r'""^o. 
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In the original notation, this imphes (37) in form of 

Vi? > 1 [ IVwPrfx < (-)"2ao f \\/u\''dx. 



{a;:|a;|>_R} 2 -'{x-:|x-|>l} 

Step 3. Rigorous derivation of the continuum version (36) for R > 1, and 
where u is now specified to be a partial derivative of the Green's function, 
i.e. u{x) = ViyG{x,x')\x'=o with i = 1,- ■ ■ ,d. In this step, as opposed to 
the previous step, we deal with the issue that we don't know a priori that 
?7^(n — u) is an admissible test function for (38). More precisely, we worry 
about the decay at |x| t oo — we don't worry about local smoothness since 
anyway, we'll apply the argument to the discrete case in the next step. As 
in Step 1, we circumvent the problem of decay through approximation by 
the periodic problem. More precisely, we consider the periodic continuum 
elliptic Green's function Gl{x,x') = Gi(a,x, x') of period L. Up to additive 
constants, it is characterized by the weak equation 

/ VC{x)-a{x)V^GL{x,x')dx = C{x')-L~'^ f C{x)dx (43) 

for all periodic Ci.^)- We note that ul{x) = Vi^x'GL{x,x')\x'=o thus is char- 
acterized by 

[ VC(x) ■ a(x)V,UL{x)dx = ViC(O). (44) 

Since VVGl distributionally converges to VVG as L '|' cxd, it is enough to 
show (44) implies 

\VuL\^dx<C{d,\)R~'^'^'' [ \VuL\^dx (45) 

f,f)^n{|x|>i?,} i[_|,|)dn{|x|>i} 

for 1 < -R < C{d)L. Indeed we can estimate the right-hand side of (45) 
using (34) and apply weak lower semi-continuity to take the limit as L — )■ oo 
on the left-hand side to obtain (36). Now, disregarding smoothness issues, 
ri^{uL — ul) is an admissible test function for (44). The argument for (45) is 
identical to the one in Step 2. 

Step 4. Rigorous derivation of (35) for R > C{d). In this step, we indicate 
the modifications in Step 2 (or rather Step 3) that are necessary to treat the 
discrete case. The first modification results from the fact that Leibniz' rule 
and thus the neat identity (39) does not hold anymore. However, we claim 
that the estimate (40) survives in form of 

d 

V{r]\u - u)) • aVn > 5^ (A(V.(r/(n - u)))^ - {[u], - u)\V^r]f^ , (46) 



22 



where [ti]j(x) = ^{u{x) + u{x + Cj)) denotes a local average of u. Indeed, 
since A < a(a;) < 1 is diagonal, it is enough to show for any two functions 
ri{x) and v{x) 

This in turn follows from the simple inequality on 4 numbers rj, fj, v, v: 

{rfv — 'if'v){v — v) — {rjv — fjvY 
= —{'q — fffvv > — (?7 — f/)^(|(f + {;))^. 

Hence, if rjlx) denotes the (slightly narrower) cut-off function for {x : \x\ > 
2R — 2} in {x : \x\ > R + 2} (which is OK for -R > 5), from (46) we obtain 
the following substitute of (41) 



a;:|a;|>2_R i=l x':_R+l<|x|<2_R-l 

< C{d,X)R-^ Yl (u-uf- (47) 

x:R<\x\<2R 

The second modification comes from the fact that we need a discrete version 
of the Poincare estimate with mean value zero on the annulus Z'^n{i? < |a;| < 
2R}, which obviously holds with a constant C{d)R? provided that R > C{d). 

Proof of Lemma 5. 

Step 1. In this step, we derive the following formulas for the continuum 

vertical derivative „ , , of spatial derivatives of the Green's function: 

^ -G{x,x') = -Vk,yG{^,y)Vk,yG{y^^')^ (48) 



dakkiy) 

d 

-Vi,,G(x,x') = -Vi,,V,,yG(x,y)Vfc,j,G(i/,a;'), (49) 



dakkiv) 

d 

V,,,V,^,'G(x,x') = -Vi,,Vk,yG{x,y)Vk,yVi:yG{y,x'). (50) 



dakkiy) 

These formulas highlight the fact that the directional spatial derivative Vi^x 
is associated with the bond b = x + cf, likewise, the continuum vertical 
derivative q^ , ) is associated with the bond e = y + Ck- In the next step, 
it will thus be more convenient to introduce the notation Vi^xC{x) = VC(&) 
and Q^ , . ( = Q^C- Equipped with this notation, we rewrite (48) - (50) in 
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the more compact form 
d 



G(x,x') = -VG(a;,e)VG(e,x'), (51) 



da{e) 

^ VG{b,x') = -VVG(6,e)VG(e,x'), (52) 



da{e) 

^ VVG{b,b') = -VVG(6,e)VVG(e,6')- (53) 



da{e 
We turn to the argument for (48) - (50): From the defining equation 

Vla{x)V^G{x,x') = 6{x - x'), 
we obtain by Leibniz' rule from applying q^ r \ that 

V>(a;) V,-^G(x, x') + V^ J{x - y) Vfc,,G(x, x') = 
oakk[y) 

and thus 

- — ^G(x, x') = -y^ G(x, z)VlJ{z - y)VkMz^ x') 
oakk[y) ^ 

= ~^b[z- y)VkM^^ z)VkMz^ ^') 

z 

= -V k,yG{x,y)W k,yG{y, x'), 

yielding the representation (48). An application of Vi^x (with i = 1,- ■ ■ ,d) 
yields (49), a further application of Vi'^ (with i' = 1, ■ ■ ■ ,d) yields (50). 

Step 2. In this step, we argue that for any two sites x and y, the depen- 
dence of the derivatives VyG{y, x) and WG{y, x) on the value a{y) of the 
conductivity is mild in the sense of 

snp\VyG{y,x)\ < exp{d^X~^) ini \VyG{y,x)\, (54) 

sup|VVG(y,x)| < exp{dh-^)mi\VVG{y,x)\. (55) 

An immediate consequence of (48) - (50), (54) and (55) is 

sup|-^G(x,x')| < C{d,\)\VyG{x,y)\\VyG{y,x% (56) 

sup|-^V.G(x,x')| < C{d,\)\VVG{x,y)\\VyG{y,x% (57) 
a{y) oa{y) 

sup|-4^VVG(x,x')| < C{d,\)\VVG{x,y)\\VVG{y,x')\. (58) 
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Coordinate-wise, (54) turns into 

d 

sup ^{Vj^yG{y,x)f 

d 

< exp{2dh~') inf Y.^V,,yG{y,x)r 

and by summation over i = 1, ■ ■ ■ ,d, (55) is a consequence of 

d 

sup ^{Vj,yVi,xG{y,x)f 
{«fefe(y)}fc=i,.-. ,d j=i 

d 

< exp{2dh~') inf 5^(V,-,V.,.G'(|/,a:))l 

{a.kk(y)}k=i,--- ,d ^~ 

In terms of the notation introduced introduced at the end of Step 1, we thus 
have to show for any pair of site x, y and any bond b 

sup V(VG(e',x))2 < exp(2t/U~^) inf V(VG(e^x))^ 
sup ^(VVG'(e',6))2 < exp(2rf5A-^) inf ^(VVG'(e', 6))^, 



{a(e)}ess " {"{e)}eGS 



e'e5 



where 5 is the set of bonds adjacent to y, i.e. {y + ei,- ■ ■ ,y + e^}. In fact, 
we shall show for any finite set S of bonds that 

max Y.{Vu{e')Y < exp{2\S\h~') mm 5^(Vu(e'))^ (59) 

e'Go e'g6 

where the field u = u{x') either denotes u{x') = G{x', x) or u{x') = VG{x', b) 
and IS"! denotes the number of bonds in the set S. At the basis of the 
argument is the a priori estimate (33) from Lemma 6, which in the edge- 
based notation reads 

^(VVG(6,e))2<A"^ (60) 

b 

We now derive the following inequality: 

^ id^) ( ^°^ ( E(^«(^'))^) ^) ) ' < sup E( VVG(e', e)r. (61) 
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Indeed, by Leibniz's rule we obtain from (52) & (53) (that in our notation 

d ' 

da{e) 

d 1 

e'G5 e'G5 

and thus by Cauchy-Schwarz 

)) 
'es e'es e'es 

which, using the chain rule, we rewrite as 



both take the form of a|^Vu(e') = -VVG(e', e)Vu(e)): 

^\ Y.^^u{e')f = ~Vu{e) J2 VVG(e', e)Vu{e') 



;-^i5^(Vn(eO)^)' < (Vn(e))^5^(VVG(e',e))2 5^(Vn(e'))^ 



^ (E(^^(^'))')')' ^ (Vu(e))2 5^(VVG(e',e))^ 



(9a(e, 

^ ' e'G5 e'eS 

We sum over e ^ S and relabel: 

E(a^(5:(v«(e')f)')^ 

< ^(Vn(eO)^sup^(VVG(e',e))2 

e'eS "^^"^ e'e5 

by the chain rule again, this can be reformulated as (61). 
Inserting (60) into (61) we obtain 



which we rewrite as 



E(5iT'os(E(V"M)f)T<A-'. 



|:(a^i°g(i:(V"M)^))'j <2A-', (62) 

We note that (62) amounts to a uniform bound of the Euclidean gradient of 
X := log ( X]e'es(^^(^'))^) with respect to the variables {a{e)}ees ^ [^, l]'"^'- 
Hence X is Lipschitz continuous on this convex set with Lipschitz constant 
2A~^. Since the set of {a{e)}ees in [0, 1]'*^' has diameter bounded by |S'|^/^, 
we obtain that the oscillation of X on this set is bounded by 2|S'|^/^A^^: 

max log ( 5^(Vn(e'))') < ^ min log ( ^(Vn(e'))') + 2\S\h-\ 

e'eS e'£S 
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which turns into (59). 

Step 3. In this step, we rephrase Lemma 6, more precisely (27), in a way 
more suitable for its application in Step 4. More specifically, we claim that 
there exists a weight exponent a{d, A) > such that 

J2\{\x-y\ + irVVG{x,y)\^'^ <Cid,X,q), (63) 

y 

for all g > 1. In fact, we claim that 

1 , , 

a := -ao (64) 

does the job. We start by noting that due to the symmetry of WG{x,y) 
under the interchange of x and y and to translational invariance (since we 
can replace a by some translation of it), it suffices to establish 

J2m + ^r^^G{y,0)\'^ <C{d,X,q). (65) 

y 

Because of g > 1, and thus f{Z'^) C f''{Z'^), we have 

J2 \i\y\ + l)"VVG(y,0)p^ <{j2\(\y\ + irVVG(y,0)p)''. 



y V 

Using a dyadic decomposition, we see 

5^|(|y| + l)"VVG(y,0)r 
y 



= |VVG(0,0)p + 5^ Yl l(|y| + irVVG(y, 

n=0 j/;2"<|j/|<2"+i 
oo 

< |VVG'(0,0)p + ^22-("+2) Y^ \VVG{y,0 

n=0 j/:2"<|i/|<2"+l 

We now may appeal to (27) to obtain 

|VVG(0,0)r + ^|(|y| + l)"VVG(y,0)|2 
y 

OO ff{A\ 

< C(d,A)(l + ^2'"("+2)2-2aon^ < C{d,X). 
n=0 
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Step 4. In this step we establish the first statement of Lemma 5, namely 
(24). More precisely, we claim that for p > max{-, 1} with a chosen in Step 
3 and all x E Z'^ it holds 



{\x\ + l)2^V(5^sup |^VVG(x,0)| 
\^^aiy) da{y) 



< C(rf,A,p)sup{(|y| + if^WVVGiyM^")]. (66) 

where C{d,\,p) denotes a generic constant that only depends on d, A, and 
p. Indeed, we first square (58) taken at x' = and sum over y: 

^ sup I -A^ VVG(x, 0) I' < Cid, A) J2 I VVG(x, y)\'\VVG{y, 0)\\ 

y a{y) Oa[Jj) ^ 

After taking the p-th power, we split the sum into its contributions over 
|y — x| < \y\ and |y — a;| > \y\ to obtain 

d „.,2\P 



(^!y.?is^^^'''^-«'i^ 



a{y) My) 

< C{d,X,p)U J2 |VVG(x,y)nvVG(y,0)|2 
y-\y-x\<\y\ 

+ ( Yl |VVG(x,y)nVVG(y,0)p)'Y (67) 

y-\y~x\>\y\ 

We first bound the first term. To this end, smuggle in a weight (|y — x| + 1)^° 
with a = a{d, A) from Step 3 and apply Holder's inequality with p and its 
dual exponent q such that ^ + ^ = 1: 

5^ |VVG(x,i/)nVVG(y,0)r)' 
y-\y~x\<\y\ 



<( E {i\y-x\ + ir\v\/Gix,y)\y') 

y-\y-x\<\y\ 

J2 [{\y-x\ + l)-"\VVG{y,0 



X 

y-\y-^\<\y 



The first term on the right hand side is bounded by Step 3, that is (63). 
After taking the expectation, we smuggle in another weight {\y\ + 1)^^*^ and 



28 



take the supremum over appropriate terms to obtain 

2p 



E ((|2/-^I + 1)"1VVG(2/,0)| 

'-A<\y\ 

( J2 (|y-^l+i)"'^"(|y|+i)"'^^)sup((|z|+i)2p^(|vvG(z,o)|2p) 



y-\y-x\<\y\ 
< 



'y-\y~'M<\y\ 
Since ly — x| < \y\ implies \y\ > ||x|, we find for tlie first r. h. s. factor tliat 

E (|2/-:,|+1)-2P«(|2/| + 1)-2W 

y-\y-^\<\y\ 

y 
Since by assumption 2pa > 2d > d, we obtain for tlie last factor 

J2{\y\ + 1)-'^" < C{d). 

y 
Combining the last three estimates yields the bound 

y-\y~^\<\y\ 

Cid,\,p){\x\ + l)-'^sup ((|z| + l)'^(|VVG(^,0)p^>^) j , 

i.e. the expectation of the first term on the right hand side of (67) is bounded 
as desired. The second term in (67) can be dealt with exactly as the first 
term by simply exchanging the roles of x and 0. 

Step 5. Like in Step 3, we rephrase Lemma 6, this time (28), in a way more 
suitable for its application in Step 6. We claim that for any integrability 
exponent q > 1 and any weight exponent /3 > we have 

sup E \i\y\ + l)"^V,G(i/, 0) l'" < C{d, A, g, /3) (68) 

a 

y 
We note that by (28) we have as soon as /3 > 0: 

j^m + ir'^^yGiyM'' 

y 



y 

oo 

<(|V,G(0, 0)1^ + ^2-''''" E \VyGiy,0)\'y'<C{d,X,l3). (69) 



y 

q (28) 



n=0 y:2"<|y|<2"+i 
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Step 6. In this step we establish the second conclusion of Lemma 5, namely 
(26). More precisely, we show that for any integrability exponent p < oo at 
least as large as in Step 3 and for any weight exponent /3 > such that 

2p{(3 -d) + d<0 (70) 

we have 

d „ ^, ^.,2\P\^ 



+ ir-'((j2^np\j—-VM^,0)\ 
^^ya(y) da{y) 



<C(rf,A,p,/3)(sup((|i/| + l)'^-i(|V,G(y,0)pP)i^ 

+ {\x\ + if-'^w^ sup ((|z| + l)'^(|VVG(;^,0)pP)^)), (71) 

for all X G Z'', where C{d,X,p, P) denotes a generic constant that only de- 
pends on d, A, p, and (3. We note that by choosing /3 small and p large, the 
exponent /3 — 1 + ^ can be made to be non-positive (in fact, as close to —1 
as we want), which proves (26). In order to establish (71), we first square 
(57) and sum over y: 

J]sup|-^V.G(a:,a;')P<C(rf,A)J]|VVG(x,?/)|2|V,G(2/,x0r 
y a(y) cia[y) ^ 

We now specify to x' = and split the sum over y: 
J]sup|-^V.G(x,0)p 
<C(d,A)( J2 + H )|VVG(x,i/)p|V,G(y,0)p 

<C(rf,A)( Yl |VVG(x,|/)nv,G(|/,0)|2 
y-\y\>^\^\ 

+ E |V,G(i/,0)nVVG(x,|/)p). (72) 

y:\y-x\>^\x\ 

We start by treating the first term on the r. h. s. of (72) in an analogous 
way to Step 4. For that purpose, let a be as in Step 3. We smuggle in the 
weight {\x — y\ + 1)" and apply Holder's inequality with p and q such that 
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P 9 



> P 
|2' 



J2 |VVG(x,y)nV,G(y,0)p 

< (5^|(|a;-y| + irVVG(x,y)p'')'"' 
y 

y-\y\>^\x\ 

The first term was bounded by a constant C((i, X,p) in Step 3. Now we take 
the expectation (■) w. r. t. a and then smuggle in a weight {\y\ + l)2p('^-i) to 
obtain as desired: 



2^' 



J2 \VVG{x,y)\'\VyG{y,0)\ 
y-\y\>h'=\ 

S/:|y|>|kl 

xsupf(|^| + l)2p(^-i)(|V,G(^,0)|2P) 



(63) 
< 



C(d,A,p)(|x| + l)-'^^'-'^sup((|z| + lfP(<^"i)(|V,G(z,0)r^)), (73) 



where we have used that 2pa > d. 

We now address the second term on the r. h. s. of (72) in a similar way, 
just exchanging the roles of VG and VVG, of y and y — x, and of a and 
—/3, where the weight exponent /3 > needs to satisfy (70). By Holder's 
inequality we obtain: 

( Yl |V,G(|/,0)nVVG(x,i/)p)' 

y-\y-x\>^\x\ 

y 
X Yl l(|y| + l)^VVG(x,y)r^ 

y-\y-x\>^\x\ 

The first term is bounded by Step 5 in form of (68). By stationarity in form 
of (|VVG(x,|/)p^) = (|VVG(x — y, 0)p^), taking expectation and smuggling 
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12^^ 



in a weight (|y — x| + 1)^^"' yields 

(( E |V,G(y,0)nVVG(x,y)p 

y-\y-A>k\^\ 
< Yl (lz/| + l)'^^(b-a;| + l)"'^''sup((|^| + l)2p'^(|VVG(2,0)|2p)). 

y:\y-x\>^\x\ 

We note that by the triangle inequality in form of \y\ < \x\ + \y — x\, in the 
range (70) the remaining sum is bounded as follows: 

J2 i\y\ + lY'^ly - x\ + 1)-'^^" 
<C{p,(3)[{\x\ + lY^^ J2 i\y - x\ + 1)-'^" 

y-\y-x\>^\x\ 
y-\y-^\>^\^\ 

Hence we have obtained 

E |V,G(y,0)nVVG(x,y)r'' 
y-\y-x\>^\^\ 
< C{d, A,p, (3){\x\ + l)2p(^-'^)+'^sup {{\z\ + 1)2P'^(|VVG(2, 0)|2P)). (74) 

z 

In view of (72), the combination of (73) and (74) and taking the 2p-th root 
yields (71). 

Proof of Theorem 1. 

We start with the proof of (3). To this purpose, we fix x G Z,*^ and p < oo; 
by Jensen's inequality, we may assume that p > po with po from Lemma 5. 
Applying Lemma 4 to ({a) = VVG'(a; x, 0) (component-wise) and inserting 
the estimate (24) of Lemma 5 yields (after redefining 6) 

(|x| + l)'^(|VVG(x,0)|2oi<C(d,A,p,p,5)(|x| + l)^(|VVG(x,0)|) 

+ (5sup (i\z\ + l)^(|VVG(z,0)|'P)i' 



We now insert (17) and take the supremum over x: 
sup ^dxl + l)'^(|VVG'(x,0)PO^) 

< C (d, X, p,p, 6) + 6 snp (i\z\ + 1)^(|VVG(2,0)|2P) 



2p\i 
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Choosing 5 = |, we obtain (3). We deal with the objection that sup^ ((|x| + 
l)'^(|VVG(x, 0)p^)^/*^^^^) may be infinite by first working with the periodic 
Green's function G^ as in the proof of Lemma 6 and then letting L f oo. 

We now turn to the proof of (4). With help of the just established (3), we 
may upgrade the result of Lemma 5, cf. (26), to 

< C(d,A,ftp)('sup ((|2| + l)"{|V,G(2,0)|*)s) + lY (75) 

Applying Lemma 4 to C = '^ xG{x, 0) (component-wise) and inserting (75) 
yields (after redefining 5) 

(|x| + l)'^-i(|V.G(x,0)PO^ 

+ ^Tsup (d^l + l)'^-i(|V,G(2,0)|2p)i) + iV 

We now insert (18) and take the supremum over x: 

sup((|a;| + l)^-i(|V,.G(a;,0)|20^) 

< C(d, A,p,p,5) + 5sup (d^l + l)'^-i(|V,G(z,0)|2p)i). 



Proof of Corollary 1. 

By stationarity, it suffices to prove the result for y = 0. It is well known 
that an LSI implies a corresponding SG, see for instance [15, Theorem 4.9]. 
Indeed, using (^^ = 1 + e/ for some /(a) in (2) and expanding to second order 
in e ^ 1 one obtains 

p ^ y <y) ^^yy) 

As in Step 2 of the proof of Lemma 4, see also [13, Lemma 11], it follows 
that 

(1/ -(/>!*> <C(p,rt((E^^^(9|y)T)- P«) 
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We fix X G Z'^ and apply this inequality to /(a) = G(a; x, 0) and use (56) 
from the proof of Lemma 5 



sup 



dG{x, 0) 



da{y) 



<C{dA)\VyG{x,y)\\VyG{yM 
a(y) ^^yy) 

to obtain 

(|G(x,0)-(G(x,0))f">'<C(rf,A,p,p)((5^|V,G(x,y)nV,G(y,0)| 

y 

The triangle inequality yields 

(|G(x,0)-(G(x,0))|'0'<C(d,A,p,p)5^(|V,G(x,y)p^|V,G(2/,0)p^ 



2\ \P 



Using Cauchy Schwarz's inequality in (■) and appealing to stationarity, we 
obtain 

(|G(a;,0)-(G(x,0)>|'^)i 
<CKA,p,p)^(|V,G(x,2/)|^^)i(|V,G(y,0)r^)^ 



2p 



2p_ 



C(d,A,p,p)^(|ViG(x-y,0)rP)i(|ViG(y,0)r^)^ 



where Vi denotes the derivative w. r. t. the first variable. Into this estimate, 
we insert the result of Theorem 1: 

(|G(x, 0) - (G(x, 0)) 1'^)' < C{d, A, p,p) Y.{\x - 2/1 + lf^'~'\\y\ + l)^^^-'^). 

y 

(77) 

We now turn to the sum on the r. h. s. of (77): By symmetry, we have 

y y-\x-y\<\y\ 

(78) 
We note that in the case oi d > 2 we have 2(1 — d) < —d so that 



J^d^l + 1)'(^~'') < C{d) < oo. 



(79) 
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Since \x — y\ < \y\ implies \y\ > \\x\ we thus have as desired for (78) 

y-\x~y\<\y\ 

y 

(79) , , 

< C(d)(|x| + l)2(i-''). (80) 

We now turn to the case oi d = 2. In this case, we spht the sum on the r. h. 
s. of (78) according to 

E = E + E 

y.\x~y\<\y\ j;:|a;-j;|<|j;| and |j;|>2|a;| j/:|x-j/|<|j/| and |j/|<2|x 

s E + E . 

y-\x~y\>^\y\ and \y\>2\x\ y:\x~y\<2\x\ and |j/|>||x| 
SO that 

Y, {\x - y\ + ir\\y\ + ir' 

y-\x-y\<\y\ 

y-\y\>2\x\ z:|z|<2|a::l 

< Ci\x\ + iy^ + Ci\x\ + iyHogi\x\+2). (81) 

Combining (80) and (81), we gather 

Yi\^ - y\ + ir'-'H\y\ + 1^'^'^ 

y 

< CM(NH-if'-^'{,„,(|J| + 2) tZl 

which we insert into (77) to obtain (5). 

Optimality of Corollary 1 for p = 1. 

In this section we will show by formal calculations that Corollary 1 is optimal 
by considering the regime 1 — A <^ 1. Recall that the Green's function satisfies 
V*a(x)VxG(x, x') = 6{x — x'). Recall that in edge-based notation introduced 
in Step 1 of the proof of Lemma 5, the components of the (diagonal) matrix 
a{x) on sites can be identified with scalars on edges via a(e) = aai^x) for 
the edge e = a; + e^. Now let a(e) = 1 + eb{e), where b is i. i. d. with 
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values at each edge taken in [—1,1]. Of course in site-based notation, b{x) 
is given by a diagonal matrix with diagonal elements b{x + Ci), i = 1, . . . ,d. 
Furthermore we assume (6(e)) = as well as (&(e)^) = 1. Note that this 
implies a G [1 - e, 1 + e] C [1/2, 3/2] (w. 1. o. g. e < 1/2), but (by linearity of 
the equation in a) all results remain true with this new upper bound on a. 
Let us expand the Green's function corresponding to a in powers of e: 

G{x, y) = Goix, y) + eGi{x, y) + . . . . 

Substituting into the defining equation for G, we find that to zeroth order in 
e, it holds 

V:V^Go{x,y) = S{x,y), 

i.e. Go is the constant-coefficient Green's function. Then to first order, it 
follows 

V:V,G'i(x, y) + V:b{x)V.Go{x, y) = 0. 

Hence it holds that 

G'i(x, y) = -Y, Go{x, z)V:h{z)V,G,{z, y) 

= - ^ V2Go(x, e)6(e) ViGo(e, y), 

e 

where the last sum is taken over all edges e. Since (6(e)) = 0, we deduce 
(Gi) = and consequently 

(Gi(a:,0)2) = ((Gi(x,0)-(Gi(x,0)))2) 

= Y. ^2^0(0:, e)V2Go(x, e')(Ke)Ke'))ViGo(e, 0) ViGo(e', 0). 

e,e' 

Since the coefficients h[x) are i. i. d. with variance 1, it follows 

((Gi(x,0) - (Gi(a;,0)))2) = 5^(V2Go(x, e))2(ViGo(e, 0))l 

e 

The behavior of the constant-coefficient Green's function is well-known, see 
for instance [16, Theorem 4.3.1], and yields that at the edge e = 2 + e,, 
the term (ViGo(e,0))^ behaves as (|2;| + l)^^'^ with a similar expression for 
(V2Go(x, e))^. Hence we find that 

((Gi(x,0) - (Gi(x,0)))2) < G{d) Y, ((k - ^1 + l)(kl + 1))'^'^'^ and 
{{G,{x, 0) - (Gi(x, 0)))^) > ttJ;^ E (d^ - ^1 + l)(l^l + 1))'^'"'^- (82) 



Gid) 
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Thus (80), (81), and (82) yield the upper bound 

If (i > 2, a lower bound can be obtained by considering only the summand 
2; = in (82). If d = 2, we restrict the sum to all z such that \z\ < \x\ to 
obtain 



((Gi(x,0)-(Gi(x,0)))2)>-^ -£ {\x-z\ + ir\\z\ + iy 

z:\z\<\x\ 



2::|2:|<|a;| 



>^(kl + l)-'log(2 + |x|). 

Thus Corollary 1 is indeed optimal in scaling. 

Proof of Corollary 2 

Step 1. Proof in dimension d > 2. First of all, the triangle inequality in 

{{■yy/''' yields 

('5^|n(x)-(n(x))|^yy''<[5^(|n(x)-(n(x))r^)'y. (83) 

^ X ' ' ^ X ' 

Since u is the decaying solution of (6) with compactly supported right hand 
side /, it can be represented via the Green's function: 

u{x) = Y,G{x,y)f{y), (84) 

V 

Consequently, an application of the triangle inequality in ((■)''p)1/(^p) yields 

\u{x) - (n(x))r^)* < Y. (|G(x,y) - (G(x, y))r^)^|/(y)|, 

y 

so that we may use Corollary 1 to the effect of 

(|n(x) - (n(x))r^)^ < C{d, \,p,r,p) 5^(|x - y\ + lY-'lfiv)]- (§5) 

y 
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We now insert (85) in (83) to obtain 

j2Hx)- {u{x))\ 

<C(ci,A,p,r,p)f5^(^^(|a:-y| + l)i-'^|/(y)iyj'. (86) 

Now let us recall the Hardy-Littlewood-Sobolev inequality in M.'^, see [17, 
Section 4.3] for a proof: 

\x-yrfiy)dyjY <Cid,a,p)njfiy)\'^dy 

for all weight exponents < a < d and for all integrability exponents 
1 < p,q < oo related byl + - = ^ + -. A discrete version can easily 
be obtained by applying the continuum version to piecewise linear functions 
on a triangulation subordinate to the lattice Z'^. We use the discrete version 
for a = d — 1, that is, 



E(E(i^-^i + i)'"'i/(^)i) J <c(rf,p)(5^i/(y)r)% (87) 

. X y / y 

in which case the relation turns as desired into - + ^ = -• Our assumption 
p > 2 and d > 2 ensure that q is indeed admissible for Hardy-Littlewood- 
Sobolev in the sense of the strict inequality g > 1. 

Step 2. Changes if li = 2. In this case, using that f{y) is supported in 
{\y\ < R}, (86) assumes the form 






■x:\x\<R 

\a;:|a;|<_R ^ y:\y\<R 

< C(d,A,p,r,p)(log^/2)f 5^('5^(|x-y| + l)i-'^|/(y)r^' 

\ X ^ y 

As in Step 1, it remains to apply the discrete Hardy-Littlewood-Sobolev 
inequality, where we note that our assumption p > 2 now ensures g > 1 even 
for d = 2. 
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Proof of Corollary 3 

Step 1. In this step, we derive the estimate 



X 

,c(,..,((E-p(EE^/te)«w)T)'- (««> 

z X y 

Indeed, it follows from the representation (84) that 
^(m(x)- {u{x)))g{x)\ 

■■5^5^(G(x,i/)-(G(x,y)))/(y)(?( 



X y 



Since the only dependence on the coefficients a is through G, the L^-version 
of SG (76), with 2p replaced by r, yields (88). 

Step 2. In this step, we estimate the right hand side of (88) as follows: 



E:s(EE^/(.M^))0 

<C(d,A,p,r)f^(^5^(|a:-z| + l)^-'^|^(x)|y(^5^(|y-z|+l)^-'^|/(2/)|) 



2\ 2 



2 ^ X 



(89) 



Indeed, expanding the square (of a vector) on the 1. h. s. of (89) and inserting 
(56) yields 

, ^ ^ dG{x,y) x2 

sup 



(EE^^/(.M^: 



-W^ . . ^«(^ 



X y 



^^ dG{x,y) dG{x',y') / ^ / /w/ w/ /x 

^ ' x,x',y,y' 

<C{d,\) J2 |V.G(x,^)||V,G(x',^)||V,G(2,y)||V,G(^,y') 

X \9ix)gix')fiy)fiy')\. 
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Consequently we obtain by the triangle inequality w. r. t. (^^ xx' yy 
(w. 1. o. g. we may assume r > 2) 



I l^/2\2/r 



(e(ee^/(.).w)7 

^ z X y ^ ' ' 

<C(d,A)| Y. (|V.G(x,^)|5|V,G(x',^)|i|V,G(^,y)|i|V,G(^,y')l^)' 

\ z,x,x' ,y,y' 



x\g{x)Mx')\\f{y)\\f{y')\ 



Holder's inequality with respect to the ensemble (■) and Theorem 1 yield 

(|V,G(x,^)|i|V,G(a;',z)|i|V,G(^,2/)|i|V,G(z,2/')l^)' 
< C{d, A, p, r) ((|x -z\ + l)(|x' - ^1 + l){\y -z\ + l){\y' -z\ + 1)) 
Hence (89) follows from partly undoing the expansion of the square: 



l-d 



( E (1= 

\ z,x,x',y,y' 



\x-z\ + iy~\\x' -z\ + iy~^{\y -z\ + 1)1-'^ 

z,x,x',y,y' 

x{\y'-z\ + iy-'\g{x)\\9{x')\\fiy)\\fiy')\ 



5^ [j^dx - .| + i)'-'\9{x)\y (j2(\y - ^1 + ^y~V{y)\) 

z ^ X y 

Step 3. Conclusion. An application of Holder's inequality w. r. t. the sum 
over z on the r. h. s. of (89) yields a bound by 

Y.[Y.^\x-z\ + IY-%{x)\^ 



2,2 

p p 
Littlewood-Sobolev inequality (87), i.e. 



with p and p such that - + i = 1 to be chosen later. We recall the Hardy- 



E(E(l^--^l + l)'"'l/(^)l)j^C'(rf,g)(^E 1/(^)1') 



1 
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if we choose p such that - = ^ + -. The Hardy-Littlewood-Sobolev inequahty 
hkewise yields 



1 1 

9 



where k = -, + k = ^ + ^ — -. Inserting these estimates into (90) and then 
into Steps 2 and 1 yields Corollary 3. 

Proof of Corollary 4. 

Step 1. We claim that for any function ri{x) supported in {|x| < R}, we 

obtain the representation formula 

{r]u){x)= ^ {u{y)VyG{y,x)-a{y)Vr]{y)-G{y,x)Vri{y)-a{y)Vu{y)). 

y-\y\<R 

(91) 

We start by noting that even on the discrete level, some aspects of Leibniz' 

rule survive, such as 

VC ■ aV(?7M) - V(?7C) ■ aVu = mVC ■ aVr] - (Vr] ■ aVu (92) 

for any function (. Indeed, since a is diagonal, (92) reduces to 

ViCVi{r]u) - Vj(?7C)ViM = uViCViT] - (Vi-qViU, 

which in turn reduces to the elementary identity 

(C - C)(^M - r]u) - {fjC - T]C){u -u)= u{C -Oiv-V)- Civ -V){u- u). 

We integrate (92): 

^ VC ■ aV{r]u) - ^ V(?7C) ■ aVu = ^ (uVC ■ aVri - CV?7 ■ aVn) (93) 
y y y 

and use it for ( = G{-,x). By definition of G, the first term on the 1. h. s. 
of (93) yields {r]u){x)] since rjG{-,x) is supported in {\y\ < i?}, the second 
term on the 1. h. s. of (93) vanishes. This yields (91). 

Step 2. We now use the representation obtained in Step 1 to obtain bounds 
on the gradient of u and consequently on the a-Holder norm of u. Specifically, 
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we claim that 

R^ sup^.|^l<^ \u[x) 



,X]\P 



^:|^l<f ?/:f<l2/l<f 

+ Rd(p-i)-P J2 Yl |V.G(y,x)r), (94) 



x:\x\<f y:§<\y\<§ 



where a and p are related hj ap = p — d. To this end, we choose a cut-off 
function r] for {|x| < -f + 1} in {|x| < -f — 1} (w. 1. o. g. R > 8). If we restrict 
to |x| < J and take the derivative in x, we obtain 

Vu{x) = J2 [u{y)VVG{y,x)a{y)VTi{y) - (Vviv) ■ a{y)Vu{y))VMy,x)), 
y 

which implies 

\Vu{x)\<C{d)R-^ J2 (|n(y)||VVG(2/,x)| + |V.G(y,x)||Vn(y)|). 

y:f<|j/l<f 

(95) 
Summing the p-th power of (95) and applying Holder's inequality, we obtain 

x:\x\<f 

<c{d,p)R-4{ Y Hy)\'Y'' E E l^^^(2/' 

y-\y\<§ x:\x\<§ y:^<\y\<§ 

+ ( E \^<y)\'Y~' E E iv.G(y,x)rY (96) 

y-\y\<§ x:\x\<f y:^<\y\<§ 

where q is the dual Holder exponent of p. Now we apply the following (dis- 
crete) Sobolev inequality: If a < 1 and p > d are related by 

a = 1 - -, (97) 

P 

then it holds that 

sup \^:^<Cid,p)( E iVn(x)r)'. 
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(This discrete version can easily be derived from its continuum version by 
extending u to a piecewise linear function on a triangulation subordinate to 
the lattice.) Therefore the left-hand side of (96) bounds the a-Holder norm 
as desired, albeit over a smaller ball. 

Let us now turn to the right-hand side of (96). It holds that 

( E l«(?/)lT <C(^,P)^'^""'^( snpjuiy)\y. (98) 

y-\y\<f a;:|x'|<f 

To estimate the second summand on the right-hand side, we note that Cac- 
cioppoli's estimate (47) implies 

J2 |Vn(y)|2<C(d,A)i?-2 ^ \u{y)\'' < C{d, \)R'^-^( sup |n(x) ' ' 
j/:|j/l<f y-\y\<R "'"'-^ 



Together with Jensen's inequality (here we need p > 2, which is obvious since 
even p > d), we obtain that 



E 



J2 iv«(i/)r)' <c{d,p)R'(^^~'^f Yl iv^(2/)i') 

y-\y\<f ^y-\y\<f 

<C{d,X,p)R^^P-^'>"p( sup \u{x)\Y. (99) 

Substituting (98) and (99) into (96) yields the claim of this step. 

Step 3. We end by bounding the Green's function to arrive at the conclusion 

sup T^^'^'-'^ 'i^': A ) < Cid, A, p, p, a) (100) 

for all a < 1, p < oo, and R < oo, where the outer supremum is taken over 
all solutions u{x) to V*aVM = in {x : \x\ < R}. Indeed, Theorem 1 applied 
to the result (94) of Step 2 yields 

^nn ^M£h:!£Mi\P 

SUPa;:|a;|<f \x\<- 

snp- '- -j—r 

<C{d,X,p,p)R"p(R'^^P~^^-P J2 J2 (1^^ - yl + 1)"^'^ 

^■■\o=\<fy:f<\y\<f 



j^ciip-i)-2p j2 Yl (1^-2/1 + 1)"^'"'^) 



x:\x\<fy:^<\y\<§ 
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if a and p are related by (97). Here we have used symmetry and stationarity 
to replace ([^/^..^(y, x)p) by (|ViG'(x — y, 0)1*') and thus bring the terms 
involving G into the form of Theorem 1. In the domains of x and y, it holds 
|x — |/| + l>|y|/2> i?/8. Therefore the first double-sum on the right-hand 
side is bounded by 

Likewise the second double-sum is bounded by 
If (97) holds, we conclude that 

\u(x)-u{Q)\ \PK 

^^P^ rT^\\ ) <C{d,X,p,p). 



Jia SUP^.|^|<^ 



uix] 



In the region {x : ^ < \x\ < R}, it obviously holds 

\u(x)-u(0)\ 8" , , ,, 

' ^ , , ^ ^' < 2—- sup u(x). 

|-^| -"' x:\x\<R 

Thus we have obtained (100) for p and a such that (97) holds. Since in (97), 
a — !■ 1 as p — 7- oo and since we can always decrease p and a in the conclusion 
(100) (in p this follows from Jensen's inequality), the estimate (100) indeed 
holds for arbitrary p < oo and a < 1 . 

Proof of Corollary 5. 

Step 1. First we write Corollary 4 in a point-wise (in a) form. The corollary 
implies that for ((■)-almost) every a, R > 0, a < 1, and all solutions u to 
V*aVu = 0, it holds that 

\u(x) — n(0)| ^,, „, 1 , / M 

sup L±J-—^Jl<K{a,R)— sup \u{x)\ 

x:\x\<R \X\ -rt. x:\x\<R 

with a constant K{a, R) such that {K{-, RY) < C{d, A, a,p) for all p < oo. 
Step 2. We argue that for every i?o > 0, there exists a sequence of increasing 
radii (-R„)„gN with lim„^oo Rn = oo such that 



N 

N 



1 ^ 

-^i^(a,i?„)<C(a)<oo. 



n=0 



Indeed, Step 1 yields that the set of random variables {K{a, -R))r>i is bounded 
in L^ ((■)). By the Banach-Saks theorem [24, Chapter 11.38], we obtain a 
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subsequence {K{a,Rn))nm with increasing, unbounded radii Rn such that 
A^""*^ J2n=o^(^^^n) converges strongly in L^((-)). Thus we have obtained a 
subsequence with the desired property. 
Step 3. In this step, we deduce the estimate 



sup \u{x) — u{0)\ < C{a) sup < -— sup \u{x)\ > 

„l-p„ ^^M lit" \x\<R„ -* 

1 ^ 1 



\x\=Ro n>M ^^n \x\<R,^ 

M 



^tl 'Rn\.\<R. 



for all A^ > M > 0. Indeed, since the radii Rn > Rq are increasing, it holds 

N 



sup \u{x) - m(0)| < ^ V sup 



n(x) — u 



Inserting Step 1 and splitting the sum yields 



sup \u{x) - u{0)\ < — Y^ K{a,Rn)— sup \u{x)\ 

1 ^^ 1 

^ ^ K \x\<Rn 

Upon noting the positivity of K{a, i?„) and the monotonicity in n of the 
suprema. Step 2 yields estimate (101). 

Step 4. We conclude by taking the limit as first iV — )■ oo and then M — t- oo. 
First, letting A^ — > oo in (101) from Step 3 yields 

sup \u{x) — u{0)\ < C{a) sup < -— sup \u{x)\ >. 

\x\=Ro n>M^Rn\x\<R„ ^ 

Since i?~" sup|2,|<^^ \u{x) \ vanishes in the limit as n — )■ oo by our boundedness 
assumption, letting M — )■ oo yields u{x) = u{0) for all |a;| < Rq. Since Rq 
was arbitrary, this implies that u is indeed constant. 

Proof of Lemma 2. 

Without loss of generality, we may assume (/) = 1. The elementary inequal- 

ity / log / - / + 1 < (/ - 1)2 then yields 

(/log/) = (/log/-/ + l)<((/-l)2). 
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Since (/ - 1)2 = {^/f- l)2(v7 + l)^ we find that 

(/log/)<((v^+lf)sup(v^-lf. 



Since (/) = 1, there exists a* G [A, l]'^ such that /(a*) < 1. Since the 
diameter of [A, l]*^ is bounded by d^^"^, it follows that 



vZ/W - 1 < a/M- vTK) < c?^ 



sup 



5v7 



da 



Likewise there exists a* G [A, 1]'^ such that /(a*) > 1 and therefore 



1 - v^Iia) < v/JM - /M <dhnp 

a 

Hence it follows that 



d^/7 



da 



(/log/) < {{./f + lf)dsnp [^y = ((v// + l)2)Zsup 
Finally it holds that 



d l/9/\2 

r^jyTa) • 



((V/ + l)2)<(2/ + 2)=4, 
and the combination of the previous two inequalities yields (16) with constant 

Proof of Lemma 3. 

The following is a simple adaptation of the usual tensorization proof, cf. [15, 
Theorem 4.4]. Take any enumeration {xn)neN of the lattice Z"^ and denote by 
(■)„, the x„-marginal of the (product) ensemble (■). Furthermore we denote 
iteratively fo '■= f and /„ := {fn-i)n- Thus /„ is the average of / over 
the first n sites. Then the 1. h. s. of (2) can be expressed as the following 
telescope sum: 

oo 

(/ log /) - (/) l0g(/) = J2 (/"-I log fn-1 - fn log fn) 

n=l 
oo 

= J]((/„_ilog/„_i)„- (/„_i)„log(/„_i)„>. (102) 

n=l 

The assumption of single-site LSI yields 

1 1 / /^f , N 2 

) • (103) 
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(/„_llog/„_l)„ - (/„-l)nlog(/„-l)n < — sup (—-^ 



Notice that the definition of /„_i immediately yields fn-i = {f)<n, where we 
have denoted shortly by (■)<n the ensemble average over the first n — 1 sites. 
By Cauchy-Schwarz, it holds that 

\da{xn)) \da{xn)'<n ^"\/ V(9a(x„)/ / <n 

Hence we find that 

(104) 

Finally we collect (102), (103), and (104) to obtain 

(/lcg/)-(/>.og(/)<l|(=upi(^)). 
which is the LSI (2). 
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